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The purpose of the present study is to show that in a particular cosmological model, with an 
affine equation of state, one can obtain, besides the background given by the scale factor, Hubble 
and deceleration parameters, a representation in terms of scalar fields and, more important, explicit 
mathematical expressions for the density contrast and the power spectrum. Although the model so 
obtained is not realistic, it reproduces features observed in some previous numerical studies and, 
therefore, it may be useful in the testing of numerical codes and as a pedagogical tool. 
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Several studies about candidates for the role of a uni- 
fied dark matter can be found in the literature flj] — [33j] - 
However, very few entirely analytical results are known, 
specially in the case of the evolution of perturbations. 
Here, the purpose is to show that, by choosing a certain 
equation of state, it is possible to have a completely an- 
alytical model, where one can obtain explicit mathemat- 
ical expressions for the density contrast 5 and, therefore, 
for the power spectrum P(k), since they are related by 
the expression [34j 



P(k) = \5(k)\' 



(1) 



Although not able to describe the current observations, 
the mathematical results shown here reproduce features 
observed in some numerical studies and, as such, may be 
useful in the testing of numerical codes and as a peda- 
gogical tool. 

The specific model studied in this text is characterized 
by an affine equation of state, where the relation between 
the pressure p and the energy density p of a cosmological 
fluid is 1 



p = po + ap , 



(2) 



defining the modified Chaplygin gas [37|-[4C 

p=( ri -l)p-Mp-», 



(4) 



where 7, M and p, are free parameters, by simply putting 
[1 = 0, M = —po and 7 = 1 + a. Since some basic 
results for the modified Chaplygin gas are available in 
the literature, it is useful to adopt such representation of 
the affine fluid. 

Using the conservation of energy in an adiabatically 
expanding universe, 



dp da 
— — = -3— 

p + p a 



(5) 



where a is the scale factor, one obtains, from eq. 



A+ (B - A) a 



-3 7 (1+a0 



(6) 



where A = Af/7 and B = pJ +A \ with po being the present 
value of the energy density of the universe. Substituting 
in this last result the value p — and using the Fried- 
mann equation for a flat space, 



where po and a are constants. Such equation is a partic- 
ular case of the quadratic model 



d\ 8tt 
a J 3 



(7) 



p = po + ap + Pp 2 



(3) 



one obtains 



with p being a constant, presented in the literature as a 
possible simple phenomenological candidate for the role 
of a unified dark matter (35j . In the affine equation of 
state, eq. ©, the constant term in the pressure can be 
seen as due to a cosmological constant, to a constant bulk 
viscous pressure or to an effective pressure associated to 
the phenomenon of particle creation [36| . This same kind 
of equation can be obtained from the equation of state 



o(t) = 



B - A\ ^ 



A 



sinh 3 ^ 



37/ftr 
2 V 3 



and, therefore, the Hubble parameter is 



H 



Sn 



A coth 



37 /8tt 



-.4 



(8) 



(9) 
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1 Throughout all the text natural units are used, i.e., c : 
k B = I- 



G = h- 



The acceleration of the universe, given by the second 
Fricdmann equation, is then 



4-7T 



(P + 3 P ) 



-in 



[(3 7 -2)p-3M] . (10) 



2 



From these results is easy to see that the deceleration 
parameter is 

a 37 1 



q ; 



aH 2 



2 cosh 2 



- 1 



(11) 



Another interesting result that one can readily obtain 
is the representation of the model in terms of a time- 
dependent homogeneous scalar field |41]. It is easy to 
show that putting 



+ , P = 



(12) 



where (p (t) is a scalar field with an associated potential 
V ((f) , one has 



In < coth 



37 
4 



8tt 



A t 



(13) 



and 



V(<p) 



— cosh 2 
2 



7 67T7 (<p - ip ) 



|2- 7 tanh 2 y/Cmj (<p - <p ) } , (14) 



where (po is a constant. 

Beyond the background given by the scale factor, Hub- 
ble and deceleration parameters, one important way of 
characterizing a cosmological model consists of studying 
the evolution of perturbations on it. Without considering 
the "averaging problem" one can write the equation 
for the perturbations (cf. equation 4.122 from Padman- 
abhan |42j|). 



15oj + 6v 2 dS k 2 v 2 8 



2a 



da 



36_ 

2^2 



H 2 a 4 

(1 - 6v 2 - 3lu 2 + 8uj) , 



(15) 



where 6 is the density contrast, u = p/p, v 2 = and k 
is the wavenumber of the Fourier mode in consideration. 
Therefore, in the family of models considered here, 

Mp -l = „.2 n .r„-l 



7-1 



Mp~ 



(16) 



i.e., v 2 is a constant, which may be negative if 7 < 1. 
Using in the equation for the perturbations the variable 



B-A 
A 



-37 



sinh 



37 
2 



8tt 



.4 



t 



and the hypothesis 

6 = x~ (I ~ x)~ 
one can obtain the equation 



(7~l)fe 2 
24tt7 2 A 



1 

37 



F(x) 
1 



(17) 
(18) 



dF 

dx 



A 



B-A 



x^^F . 



(19) 



A solution of the equation ([19]) is easily found only for 
a few specific cases: k = 0, 7 = 1 and 7 = 2/3. The 
last one is more interesting because then one obtains a 
density contrast which depends on the wavenumber k. 
Explicitly, the solution for this case is 



S = aSi + c 2 <5 2 , 
with ci and c 2 being arbitrary constants, 



(20) 



51 = 
and 

5 2 = 



x 2 



(i- x y 



■ 2-^1 



(21) 



— 2F1 



k 2 



13 
4'2 



;x 



(22) 

where A; 2 = k 2 / [32tt (B — A)]. Another equivalent way 
of writing such result is 



5 = 



x 2 



(v arcsin x 2 ^ + c 2 sin (v arcsin x 2 J 



where z/ 



ci cos ( v arcsin x 
1/2 

1 + 4fc y ' 



(23) 



There is yet a third way of writing the solution given 
in equations (f2"0"]) , ([2"Tj) and (|2"2"1) , which is obtained from 
the aplication on equation ([23]) of the definitions of the 
inverse trigonometric function arcsine L 43| and of the vari- 
able x. Then one has 



c' 2 S- 



S = c[S + 

where c[ and c 2 are arbitrary constants, 
1 / cosh i =F 1 



S± 



and where v 



sinh t cosh 2 t \ sinhi 



(24) 



(25) 



1 



.1/2 
4fc 2 ] and t 



(47rM) 1/2 t. This 

result can be obtained directly from the equation for the 
density contrast using t as variable [42| . 

v 2 

5 + H[2-3(2uj- v 2 )} S + k 2 —S 
3 



-H 2 (l-6v 2 -3uj 2 + 8lj) 8. 



(26) 



Using eqs. ©, ©, ©_and ([16]), with 7 = 2/3, changing 
the variable from t to t, and choosing 



S = 



G(t) 



sinh t cosh 2 t ' 



one obtains the equation 

d 2 G , -dG 
coth i 



rfi 2 ' " d* sinh 2 i G ' 



(27) 



(28) 
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FIG. 1: Curves showing the variation in time of the function 
5+ (fc) for some values of k. All curves present peaks whose 
location and magnitude varies with k. Notice that the time 
variable is t — (4nM)^ t, where t is the usual cosmological 
time, with M being a parameter of the model (see eq. Q). 
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FIG. 2: Curves showing the variation in time of the function 
8- (fc) for some values of k. Notice that the time variable is 

t = (4-7rAf ) 5 t, where t is the usual cosmological time, with 
M being a parameter of the model (see eq. (|4|l). 



which admits as solutions the toroidal functions 
Pq (cosh t) and Qq (sinh t) [43| or, simply, the result given 
by equations ([24]) and (pS| . 

Figures [T] and [2] show the behaviour in time of the 
functions S + and <5_ for some values of k. It is easy 
to notice that 5- represents rapidly decreasing modes, 
which dominate for smaller times, while 5 + represents 
modes with peaks which become the dominating ones 
for greater times. The presence of peaks in the density 
contrast is a feature shown in previous numerical studies 
@; EH an d' therefore, confirmed here analitically. 

It is easy now to construct the power spectrum P (fc, i) 
associated to this solution, with two types of modes. Es- 
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FIG. 3: Graph showing the behaviour of the power spectrum 
P, as given by equation (|29|) . for some values of k, at the 
time t = 0.3, chosen arbitrarily. Notice that the scale in the 
vertical axis is arbitrary. 
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FIG. 4: Graph showing the behaviour of the power spectrum 
P, as given by equation (|30[1 . for some values of k, at the time 
i = 5, chosen arbitrarily. Notice that the scale in the vertical 
axis is arbitrary. 



sentially, if the decaying mode is dominating one has 

p(M)~^[,M0] 2 , (29) 

while if the growing mode is dominant one has 

P(k,t)^cl[5 + {t)f . (30) 

Figures [3] and |4] show the power spectrum thus obtained 
for arbitrary times. 

It must be clear that none of the results shown here 
are consistent with the current observations. For exam- 
ple, they are not compatible, at the same time, with an 
initial scale invariant power spectrum (P oc k) and the 
effectively measured power spectrum [44[ : if one sets the 
constants c\ and C2 (or d x and c' 2 ) to give the correct 
initial conditions, one misses the target of real data, and 
vice-versa. Therefore, the model shown here is simply 



4 



a toy model, or a mathematical possibility, but an in- Acknowledgements 
teresting one because it has the peculiarity of yielding 
analytical results along all the way. Anyway, this may be 
seen as a practical example of the idea behind the famous 
Einstein's phrase "As far as the laws of mathematics re- 
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